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The article deals with the problem of finding integrable boundary conditions for multi- 
dimensional equations. The research is stimulated by the following problem, posed by Richard 
Ward: Find discrete versions of the well known generalized Toda chains, corresponding to the 
Lie algebras of finite growth (see R.Ward, [I]). Really, one has to find cutting off constraint for 
the infinite discrete Toda chain 
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Q \ where / = f{u,v,k) - is the field variable, and sub- and super-indices denotes shifts in the 

^ ! corresponding arguments: = f{u + 1, f , k), fy = f{u, f + 1, fc), = f{u, v,k + 1), /"'^ = 

§: f{u,v,k-l). 
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But in the paper we study the other well known model - the Hirota equation 

TuT, - TTu. = T^T-\ (2) 

which is also a discrete version of the Toda chain. However the algorithm proposed can be 
applied to any kind (discrete and continuous [?) of integrable equations with 3 independent 
variables, including the Ward's problem. Cutting off constraint for Hirota's equation are found 
corresponding to the generalized Toda chains of the series and C„. 
^ ■ We tried also to answer the question: What is the Lax pair of the boundary value problem? 

CN ' It differs from that of the equation. Our observation is that boundary condition put on the field 
variables generates some boundary conditions for the eigenfunctions (solutions to the linear Lax 
^ \ system). As a rule one has to take two different Lax pairs with the eigenfunctions related to 
O \ each other at the border points. The best illustration to this idea provides the example of 
, the Kadomtsev-Petviashvili (KP) equation below in the section 4. It can be checked by direct 
^ ■ computation that the boundary problem ((23), (IHUj) admits the Lax pair (jHTj) . (jHH). In 

^ ■ other words the boundary problem ()35|). ()36p is a consistency condition of the linear boundary 
^ ' value problem (j37jl . (j38|) with the time evolution (j39|l . It is remarkable that passing in the 
\ aforementioned example from the eigenfunctions to the Gel'fand-Levitan-Marchenko kernel 
leads to a linear boundary value problem which is solved by separation of variables. 

1 Algorithm of searching for BC's 

How to extract the integrable BC's by using only the Lax representation of the equation? 
To answer the question let us concentrate on the example of the Toda chain 

Uxt{n) = exp{u {n — 1) — u (n)} — exp{u (n) — u{n + 1)}, (3) 
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lis L,ax pair is aennea oy ijapiace xransiormaLions lor ine loiiowmg nyperDoiic equation 



0xi (n) = -u, {n) 0, {n) - e«(«-i)-«W0 (n) , (4) 

ct){n + l) = {D., + u.,{n))<P{n), (5) 

0,(n) = -e"('^-i)-"("V(n - 1). (6) 

Take one more Lax pair of the same chain essentially different from the above one 

V-xt {n) = -ut (n) (n) - e"^""^)-"^"^ (n) , (7) 

^(n + 1) = (A + nt(n))^(n), (8) 

^^(n) = -e"("-i)-"(")^/>(n - 1). (9) 



Proposition 1 (see j2]). Let the operator of the form M = aD^ + bD'^ + cD^ + d exist such 
that for n = for each solution ip of the equation Q the function = M-^ solves the equation 
(^. Then the field variables m(0), «(—!), m(— 2) are connected by one of the constraint 

1) e"(-i) = 0, 

2) «(-!) = 0, 

3) «(-!) = -n(0), (10) 

4) u^{-l) = -w,(0)e-"(°)-"(-i), 

5) e"(-2)-"(-i)(n.(-2)+M,.(-l)) = 
= e«(-i)-«(o)(6w.(-l)-2M.(0)). 

The corresponding differential operators M are of the form 

1) M, = aoe''{Dl + 2u,Dl + {u^, + ul)D,) + 
+boe''{Dl + u,D,) + coe'^D,, 

2) M2 = e"D2 ^ M,e"D,, 

3) M3 = e"D,, (11) 

4) M4 = e^D^ + M^e"D^ + e"", 

5) Ms = e"D^ + 2u,e"D^. + 

+ (^ixx - u^x{-l) ul{-l))e''D^. 

where a^, bo, Cq - are real constant parameters, and u = u{0). 

Conclusion: Under integrable boundary conditions two essentially different Lax pairs (which 
are not related to each other by gauge transform) become conjugate along the border. The 
conjugation relation is nothing else but the corresponding boundary condition for the eigenfu- 
nions. 



2 



^ JJUUllUdl^ I^U11U.1LHJ11» lUl Lilt; XVJT tJl^llCltHJll I^Ull»l»tt;ilL 

with the Lax pair 

Apply the method suggested above to the KP equation 

Ut+Uxxx-QuUx = -So^Wy, (12) 

Wx = Uy. 

Remind the Lax pair for it which is defined as the following system of the linear equations 



= a(j)y + u(j)^ (13) 
-^xxx + Qu(j)x + 3(ux + aw)(j). (14) 



The dual pair is of the form: 



4>xx = -a(l)y + u4), (15) 
04 = -4:(t)^xx + Qu4>x + 3{ux - aw)(j). (16) 

Lemma 1. Solution u = u{x,y,t) to the KP equation satisfies the boundary condition 
w\y=o = (or the boundary condition Ux + aw\y=o = 0) iff for each solution = (j){x,y,t) of 
the equation (fT^ the function = (p\y=o (respectively the function = Dx4>\y=o) solves the 
equation (fT^ for y = 0. 

3 Zakharov-Shabat dressing and the boundary problem 

Let Uo = be a starting solution to the KP and let 0o = (poi^y t) be a starting eigenfunction. 
Find a new (dressed) eigenfunction = (j){x,y,t), corresponding to the new (dressed) potential 
u. To this end we use the following integral operator with Volterra kernel K = K{x, z, y, t) (see 

my- 

(t){.x,y,t) = 0o(a;,?/,t) + 

+ r K{x,z',y,t)Mz',y,t)dz' (17) 



Its kernel K satisfies the Gel'fand-Levitan-Marchenko equation 

K{x,z,y,t) + F{x,z,y,t) + 

+ r K{x,z',y,t)F{z',z,y,t)dz' = 0. (18) 

J — OD 

The kernel of the Gel'fand-Levitan-Marchenko equation F satisfies in turn a system of linear 
differential equations with constant coefficients 

^dF^_d^^^^Q ^^g^ 

dy dx'^ dz'^ ' 
dF ^^d^F d^F^ , , 
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iiie searcnea soimion u oi lue nonimear equation is aennea xnrougn xne Kernel n oy means oi 
the formula 

u{x,y,t) = 2—K{x,x,y,t). (21) 

Formulate now the dressing method based on the other Lax pair, then the system p9|) - ()20|) is 
replaced by the following system 

_ ^ ^ ^ ^ Q ^22) 

and the kernel K of the transformation operator, which converts the starting eigenfunction 0o 
into an eigenfunction of the new associated equation 



4){x,y,t) = (f)o{x,y,t) + K{x,z',y,t)(f)o{z',y,t)dz' (24) 

J — oo 

will solve the similar Gel'fand-Levitan-Marchenko equation 

k{x,z,y,t) + F{x,z,y,t) + 

+ r K{x,z',y,t)F{z',z,y,t)dz' = 0. (25) 

J — oo 

Theorem (see j2j). Involution 

dF OF 

— {x,z,0,t) = - — {x,z,0,t) (26) 
corresponds to the boundary condition 

+ aw\y=o = 0, (27) 

and the involution 

F{x,z,0,t) = F{x,z,0,t), (28) 
corresponds, respectively to the boundary condition 

w\y=o = 0. (29) 

Consider in more details the KP2 equation, putting in all formulae above a = 1. Let the 
Gel'fand-Levitan-Marchenko equation have the degenerate kernel 

N 

F{x,z,y,t) = J2^new{{ql-Pl)y -4:{(ll+pl)t + qnZ + Pnx). (30) 



Then evidently 



n=l 



N 



F{x,z,y,t) = J2cnew{{(ll-Pl)y -4:{(l'n+pl)t + qnX + Pnz). (31) 



ra=l 
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ouppose 111 aaaiLion inax lue Kernels r ana r are reiaiea oy ine mvoiuxion ijzpji . wnicn can oe 
represented as 



N 



CnPn exp{g„2; + p„x - 4:{ql + pl)t} = 

N 

CnPn exp{g„x + PnZ - 4(g^ + pl)t}. 



n=l 

N 



n=l 

It implies immediately that 

Pn = QN-n, CnPn = -CN-nQn- (32) 

Thus, under constraint ()32|) the solution of the KP equation given by the formulae 

u{x,y,t) = -2— In det A, (33) 
where A is a matrix with the entries 

A — X _L. J(q?,-pi)y-MQ^+P^)t+(Pn+qm)x) (n4\ 

Pn + qm 

satisfies the boundary condition (j^7j) . which is + w\y=Q = for the KP2. Its differential 
consequence with respect to x looks more illustrative: Uxx + Uy\y=o = 0. 



4 Boundary problem on the stripe < y < 1 

Consider the boundary problem on the stripe < y < 1, 

Ut + Uxxx - Quux = -So^Wy, (35) 

w\y=o = 0, w\y=i = 0, (36) 

which is equivalent to a compounded Lax pair. Eigenfunction solves the boundary problem 

with boundary conditions imposed along the same direct lines y = and y = I 

(l)-^\y=o = 0, <f>-^\y=i = 0. (38) 
Due to the Lemma 1 the conditions ()38|1 are consistent with the time evolution 
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Accoramg lo tne tneorem aoove tne iviarcnenKo Kernel wnicn is not; scalar now oui a vecxor- 
function satisfies the boundary value problem 

-^y{p)=^'^{^p)-^'^{^p) (40) 

with the conditions 

F-F\y=o = 0, F-F|,=i=0. (41) 
Concentrate on the case of the degenerate kernel 

N 
n=l 

By means of the constraint F{x, z, y, t) = F{z, x, y, t) one gets 

N 

F{x,z,y,t) = J2(^nexp{{ql-pl)y -4:{ql+pl)t + qnX + PnZ). (43) 

n=l 

Substituting the boundary conditions in the points y = and y = I one gets 

C„ = CN-n, Qn = PN-n, PN-n = Pn, pI - pI = '^ikn. (44) 

The last constraint shows that the quantity p„ is described by an entire and a real parameters 
G Z and /3„ G R such that 

/ , N k^ 

Pn =Pn{k) = ITTT+Vn 

for n = l,...N/2. 

One can replace each of the conditions (jHSl) (or both simultaneously) by the condition 
Ux + w\y=yo = taking either y^ = or yo = 1. Then the boundary conditions for the Lax pair 
(^Hjl will be replaced by — ipx\y=o = 0, or respectively, (p — ipx\y=i = 0. The Gel'fand-Levitan- 
Marchenko kernel will satisfy instead of (jlT|l the following boundary condition Fx + Fz\y=o = 0, 
or respectively, Fx + Fz\y=i = 0. 

Remark. The problem is open how to apply the inverse scattering transform method to 
the problems on the stripe. 

5 Cutting off constraint for discrete chains 

The well known Hirota equation is a reduction of the following basic equation (see, for instance. 

El) 

(A-l)(^^^^J-^) = 0, (45) 

Here t = t{u, v, k) is the dependent variable. Super- and sub-indices denote shifts with respect 
to discrete arguments: t~'' = t{u, v,k — l),t^ = t{u, v, k + 1), t„ = t{u+l, v, k),tjj = t{u, f + 1, k). 
In order to get Hirota equation one should put the expression inside of the curly brackets in 
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i|40|i equal lo —i. ine L,ax pair lor tne equation ij40ji is given as pair oi cnams oi ijapiace 
transformations 

^u = ^' + (46) 

= ^ + (47) 
for the following discrete hyperbolic equation 

i^uv = ^u + ji^A - If ^, (48) 

where ipu = 'ipiu + fc) and ip^ = ip{u^v + l,fc) etc. (about discrete Laplace transforms 
see survey by I.A.Dynnikov, S.P.Novikov, |5J. The iterations of Laplace transform above are 
enumerated by upper index. The operators forming these transforms are Du and D^. There is 
one more Lax pair for the equation ()45j] . which is based on the other kind of the Laplace chain 
generated by operators and 

tt^ 



Lk 



+ (49) 



lk +—k 



serving the folloving hyperbolic equation 



4. 4.k 4. 4. 4. 4.k 

(Puv = 7n—(Pu + 11— (Pv - T^(P- (51) 

Note that discrete hyperbolic equations admit two more pairs of mutually inverse Laplace 
transforms based on the operators D^^, Dy and Du-, D^^. So totally they have four kinds of 
Laplace transforms. Two pairs define as one could see above two different Lax pairs. We show 
below that two others define the cutting off constraint. 

Proposition 2. Let the equations pHjl and (jHTjl be connected by substitutions 

i) = {aDu + 6)0, 

then the following boundary condition (cutting off constraint) is satisfied 

And the substitution and its inversion are of the form 

V' = |(l-/^«)0, = ^(D;1-1)^. (52) 

These substitutions are nothing else but the mutually inverse Laplace transforms with op- 
erators Du and -D~^, respectively. By replacing u ^ v one would get the other kind of BC's. 
Let us reduce the chain into a segment k G (0, A^) by setting the BC's at the ends /c = 0, 
k = N: 

t(M,w,-l) = 0; {Du - 1) """^^fc^a" " lfc=^ = Q' 
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mgenmncxions (p[u,v,K) ana ip[u,v,K) are aiso restricxea lo me same segmem k fc (,u,iv j. j\i 
the left end the hnear equations of the Lax pair are cut evidently under BC. At the right end 
one has to prescribe any of two equivalent equations (jK^ taken at the point k = N. 

Proposition 3. Let the equations (jIHj) and be connected by second order substitutions 

tjj = (D„ - - 1)0, 

=iD;'-l)ziD;'-l)i;, (53) 

defined as compositions of two Laplace transforms, where z = then the following boundary 
condition is satisfied 

The substitutions found allows one to cut up the infinite Lax pairs and adopt them to 
restricted chains. 

6 Cutting off constraint for Hirota's equation 

Let us pass to the Hirota equation 

T,T, - TTu. = T^T-\ (54) 

For this reduction the boundary conditions above can be integrated under which they take (up 
to the point symmetries) the following forms 

i) T(u,t; + l,l) =T(M + l,t;,-l); 

ii) T(u,t;,0) = T(u + 1, w, -1). 

In the continual limit the Hirota equation approaches the Toda chain (see, for instance, |3]), 
the boundary conditions found above turn into cutting off constraint for the chain corresponding 
to the Lie algebras of series C„ and Bn- The case corresponding to the series An is well known, 
it has been found years ago. But the case of Dn is not found yet, it is connected with rather 
labour consuming computations. 

The discrete Toda chain (P) can also be studied by the method above. 
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